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Abstract 

The low energy photoemission spectra in quantum antiferromagnets are stud- 
ied by using several approximation- free calculations and rigorous theorems. 
The important and measurable property found is that the hole eigenstates 
with momenta differing by the antiferromagnetic wavevector Q are equivalent 
and degenerate in energy. However the corresponding eigenstates differs by 
the presence or the absence of a well defined quasiparticle corresponding to 
a singular -zero energy- magnon, carrying spin one and momentum Q. This 
difference between the two eigenstates affects dramatically the spectral weight 
as a function of the scattered momentum, since a sharp effect at the surface of 
the magnetic Brillouin zone is predicted, in apparent agreement with recent 
experimental data. 
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I. INTRODUCTION 



Soon after the discovery of high-T^ materials it was proposed that, this kind of interesting 
phenomenon , could be explained by the strong antiferromagnetic fluctuations characterizing 
the phase at small hole doping. The fact that in this regime the materials are close to a 
Mott-Hubbard transition have stimulated an intense scientific activity on the study of simple 
one band models , like the Hubbard model and the t — J model, whose apparent simplicity 
is very attractive. However these models represent a still unsolved many-body problem in 
condensed matter physics, as it is still questioned whether they can have a superconducting 
ground state in two spatial dimensions. 

The so called t-J model is defined by the following hamiltonian: 

H = -t i4aCja + h.c.) + JY1 (S' ■ Si - ^niUj). (1) 

<i,j>,a <i,j> 

where the constraint of no double occupancy is understood. Here c|(cj) creates (annihilates) 
an electron at site i, rii = J^^ia is the corresponding density operator with riia = cj^Cia, 
the symbol < i,j > means summation over nearest neighbors, J is the super-exchange 
coupling and finally the spin density operator Sj is defined by the Pauli matrices a: Sj = 
J2 4,cr~^Cj^o-'. Hereafter periodic boundary conditions are assumed on a lattice containing 
L sites with holes {L — Nh electrons). In these paper I will mainly focus on the property 
of a single hole Nh = 1, which was a subject of a vast literature starting from Nagaoka 
and Brinkmann and Rice in the early seventies. This problem is also important because 
the property of a single hole in a quantum antiferromagnet are experimentally accessible by 
angle resolved photoemission experiments that quite recently have reached an high level 
of resolution and reproducibility especially in two dimensions. In fact , in absence of hole 
doping, it is widely accepted that the Copper-oxide planes are well described by an effective 
Heisenberg model ( the term proportional to J in |1|) with long range antiferromagnetic order 
in two dimensions . When an electron is photoemitted , a single hole is free to propagate in 
the Copper-Oxide plane, by hopping from Copper site to Copper site with matrix elements 
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described by the t term in Eq.(||). It is reasonable to assume that interaction with the p 
Oxygen bands and other effects cannot lead to a qualitative change of the model provided 
the Coulomb repulsion remains large to forbid double occupied sites and properly define the 
position of the hole in the lattice. In fact by taking the large U limit of the Hubbard model, 
the lower Hubbard band is correctly described by the t — J model, which represents therefore 
a realistic model for the interpretation of photoemission experiments at low energies. 

One of the most important property is to understand whether a single hole can propagate 
as a free particle as in a conventional band insulator. In this case if an electron is extracted 
from the lower band it will remains for infinite time with the given hole momentum k. 
Within the quasiparticle picture it is clear that interaction can provide a finite lifetime, but 
since in an insulator all the possible states in the lower band are occupied there is no phase 
space for decaying processes at least for the lowest possible excitations. As in the Landau 
theory for Fermi liquids it is then possible to define the so called quasiparticle weight Zp 
which measures the distance between the many body state |p > of a single hole and the 
corresponding ground state of the insulator \H > in which a hole has been created with 
momentum p and spin a, Cp^o-\H >. The quasiparticle weight is then given by: 

Zp = \ < p\cp,^\H > \^ (2) 

For a band insulator with no interaction Zp is one by definition, whereas if a quasiparticle 
picture is well defined Zp remains finite in the infinite volume limit. A key question is to 
understand whether in two spatial dimension a quasiparticle theory can be defined for an 
insulator or if this quantity vanishes with a power law as in one dimension In particular 
in an antiferromagnet it is crucial to understand if the scattering with the underlying low 
energy modes-the spin wave excitations- can lead to some anomalous effect. 

By the Lehman decomposition of the one particle Greens function Gp{u!) the finiteness 
of the quasiparticle weight Zp determines a delta— function in the spectral weight A{uj < 
0) = —^Gp{uj), which , as mentioned, is experimentally accessible by ARPES. However 
it is difficult to distinguish , in an experiment, which is biased by the resolution of the 

3 



electronic device, if a true S function exists or if the spectral weight is completely incoherent 
, i.e. without 6 function contributions. This 6— function weight Zp occurs at an energy Ep, 
which represents , within the quasiparticle picture, the lowest one hole energy with given 
momentum p, and Ep properly describes the energy dispersion of the hole. 

In two spatial dimensions there is no exact solution of the "one hole problem" but there 
is a general believe that the basic dynamical properties of a single hole in a quantum anti- 
ferromagnet are well described in the paper by Kane Lee and Read (KLR) f^/,. In this work 
the problem is solved using two main approximations: a large spin limit and a self consis- 
tent Born approximation to sum the relevant diagrams contributing to the Greens function, 
accounting for the scattering of the hole with the spin-waves. Within these approximations 
i) the quasiparticle weight is finite , ii) the minimum possible energy Ep is generally at 
momentum p = (f, |) and finally , as a general property coming from the long range an- 
tiferromagnetic order, iii) momenta differing by the wavevector Q = (vr, vr) are equivalent 
both for the energy and for the quasiparticle weight: 

Zp = Zp^Q 

Ep = Ep+Q. (3) 

In the original KLR paper there was some speculation that Zp could vanish away from the 
bottom of the band even within their approximations. However it was shown later that 
Zp is strictly finite for all possible momenta in the Brillouin zone as long as J > 0. 

An accurate determination of the single hole dispersion was recently determined by 
Quantum Monte Carlo and is essentially in agreement with the KLR result. The minimum 
of the band at p = (f , |) is also confirmed experimentally. For the quasiparticle weight 
Quantum Monte Carlo is not accurate enough and , quite recently, the accuracy of the KLR 
theory has been questioned since numerical data obtained with Lanczos exact diagonalization 
up to 26 sites clearly showed that Zp^q « Zp ai least for p ~ 0. 

I will show in this paper (and prove for J = 2t) that while the spectrum satisfies ( 
rigorously shown in App. [VIII| using a variational ansaz), the quasiparticle weight Zp ^ Zp^q. 
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P] This effect induces a Icind of Fermi volume discontinuity which is consistent with recent 
experiments by Wells. In this work it was pointed out that the photoemission experiments 
of insulating antiferromagnetic materials look similar to the corresponding metallic ones. |1| 
This simple effect is instead a direct consequence of a singular magnon excitation present 
in a quantum antiferromagnet which cost zero energy and change the spin by one and the 
momentum by Q. In essence this strange and measurable effect is intimately related to the 
non commutativity of the order parameter m with the total spin S"^. In any mean field 
theory one usually selects a direction of the order parameter and in this way the spin is no 
longer defined for the approximate ground state. In an experiments however the variation of 
the total spin is measured since the outcoming electrons have a definite spin. An approach 
which is capable to work in a basis with definite spin S"^ and order parameter m? allows 
to deal correctly with the excitations and to determine the sharpness of the photoemission 
spectra in a quantum antiferromagnet. 



II. FORMALISM 

I consider the t — J Hamiltonian (0) and follow the formalism recently proposed by 
Angelucci et ai. f^. The basic advantage of this approach is the possibility to simplify the 
local constraint of no doubly occupied sites by extending the hamiltonian to an enlarged 
space, where the physical subspace, satisfying the constraint, is obtained when the eigenvalue 
of a spin-operator L is maximum. In this representation the t — J Hamiltonian matrix reads 



H = Eo+Yl t Ui+rjR + ^(1 - - nn^^^ 



Xr,r+t.- (4) 



where /j^ and fR are canonical spinless fermion operators creating or annihilating a hole at 
site R and commuting with the corresponding spin operators Sr. Moreover ur = fjifR is 
the hole number at the site R, Nh = being the total number of holes, and the operator 

R 

Xr,r' = 2Sr Sr' + ^ (5) 



interchanges two S = 1/2-spins at sites R and R'. Hereafter periodic boundary conditions 
over a finite hyper-cubic box with hnear size I and number of sites L = l'^ is assumed. The 
constant shift of the energies Eq = ^{2Nh — L), which represents the classical contribution 
of an antiferromagnetic Neel state, is also omitted. 

The hamiltonian is naturally defined in a space where there is a single spin for each site 
and the hole can be created over each spin site by However the spin value at the position 
of the holes is clearly unphysical because in the physical Hilbert space the hole does not 
carry spin. As anticipated this difficulty can be easily solved because the hamiltonian leaves 
invariant the physical subspace of states where the spins over the sites occupied by the holes 
are all frozen to 1/2 (or —1/2) . In fact the following pseudo-spin operator, measuring the 
total spin at the positions of the holes, and satisfying the algebra of the angular momentum 

L = ^tirSr, (6) 

R 

commutes with the Hamiltonian and the invariant subspace with maximum spin L = 
Lz = coincides with the physical Hilbert space of the t — J model, where the two models 
have the same matrix elements , i.e. they identically coincide. 

Analogously the physical total spin has to be measured where there are no holes and 
within the present formalism it is defined by: 

S = J2il-nn)Sn, (7) 

R. 

Even this operator commutes with the hamiltonian and obviously with the pseudo-spin L. 
In the following sections the total spin 

S'°' = S + L (8) 

will be introduced too. 

Any eigenstate |^ > of the hamiltonian can be labeled by six quantum numbers: 

• The number of holes Nh 
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• The value of the physical spin S and its component on the z-axis Sz 

• The value of the pseudo-spin L and its component on the z-axis 

• The total lattice momentum p of the state. 

and the physical sector is obtained when L — — Correlation functions are expectation 
value of operators O over these eigenstates < ilj\0\ip > . If such an operator commutes with 

— * — * 

the total spin S or the pseudo-spin L -as, say, the hamiltonian- the corresponding expectation 
value does not depend on the spin component on the z-axis Sz or L^. In particular all the 
eigenstates differing by the quantum number Lz or Sz are degenerate. 

At the end of this section we remind that a state \ipp > with definite lattice momentum 
p is an eigenstate of the translation operator Tr, such that Tji\t/jp >— e^'^^lV'p > The 
operator Tr is the translation operator that brings the origin O to the lattice point R. It is 
formally defined by the following relations valid for any R and R': 

TRfR'T-R = Jr+r' 

TrSr>T_r — Sr+RI (9) 

The lattice momentum p is defined in the Brillouin zone (BZ) with each component \pi\ < tt 
(the lattice constant is one). An important region contained in the previous one is the 
magnetic Brillouin zone (BZ') defined as the locus of points such that 7p > where jp = 

III. TOWER OF STATES IN A QUANTUM ANTIFERROMAGNET 

In a quantum antiferromagnet the most important quantities for the correct description 
of the low energy physics are the order parameter rh and the total spin S, that, within our 
formalism, read: 

m^We^'^^SRil-riR) 

^ R 

S = J2SR{l-nR) (10) 

R 
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where Q is the antiferromagnetic wavevector Q = (vr, vr, ■ ■ ■) 

The commutation rules of these two operators are well known 0: 

[Sj, Sk\ = i€j^k,iSi 
[mj,mfc] = j^ej,k,iSi (11) 

In particular using the above relations it is easy to obtain the commutation of rri^ = 
rUx + imy or = irix — imy and with the total spin S"^: 



-2m,5+ + 2m+(l + 5,) 



S'^.m =2mzS + 2m {1 — 



S 



m S'^ — m^S + 2m, 



(12) 



The operator commutes with S"^ and Sz- In the following it will be assumed that 
also commutes with the hamiltonian since if long range order sets in becomes a 



macroscopic classical variable in the infinite volume limit. |]IU| Then we may classify the 
eigenstates in terms of the four quantum numbers S*^, 5*2, and the total lattice momentum 
of the eigenstate. 

Suppose that an eigenstate \S^rn?,f) > is given with total momentum definite spin S 
and maximum azimuthal spin Sz = S , (all the other 5*^ component can be easily obtained 
by a repeated application of the total spin lowering operator S~), it will be shown in the 
following that the state m^\S, rn? > has spin S' + l, Sz = + the same order parameter 

but momentum p + Q. 

The change of momentum is easily understood since TRiriT^R = e^^^^m, yielding 
Tjim~^\S,m'^ ,p >= e'^^~^'^^^m'^\S,m^,p >, i.e. a state with momentum p + Q. On the 
other hand, by hypothesis: 



S^\S,m\p>= S{S + l)\S,m,p> 
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whereas by applying (|T2D: 

S^m+\S,m'^,p>= {m+S^+ [S^,m+]) \S,m^,p>= + 1)(5 + 2)m+|5, m^, p >, 

i.e. m'^\S,rn? ,p > is a state with definite spin equal to 5* + 1. Analogously Sz is defined 
since, by [S'^,m"'"] = m"*", it follows Szm'^\S,rn? >= {S + l)m^\S,m'^,p >. Finally 
does not change by applying m"*" to the state jS", m^,p > as [m^,m+] = 0. 

By the simple assumption that commutes with the hamiltonian for L oo, it has 
been derived that all the different spin sectors are degenerate and the momentum is defined 
modulo Q. In a finite system this degeneracy is slightly removed and in fact the energy 
spectrum as a function of the spin behaves as a free quantum rotator spectrum determined 



by the value of its static susceptibility xst jlTiP^ 



Instead by spin-wave theory the magnon excitations are characterized by much larger energy 
costs ~ y. It seems therefore exact at least at low enough energies that the tower of states 
with different spin can be considered degenerate for L oo, when becomes classical. In 
a bipartite lattice without holes the spin is integer and the lowest eigenstates with S = n 
can be obtained by applying n times the operator to the singlet spectrum of eigenstates 
\S = 0, m^,p >: 

\S + n,m'^ ,p + nQ >(x {'m~^)"'\S,m'^ ,p > (13) 

In particular, from 2Q = modulo 2it, odd integer spin have the lowest energy state 
with momentum Q, whereas even spin have vanishing momentum (referenced to the singlet 
momentum which is or Q according to the parity of ^ respectively). 

In presence of a single hole or an odd number of holes the spin is half odd integer (for 
half odd integer spin s) and the minimum one is S* = |. Consequently the eigenstate with 
minimum spin, generator of all the tower of states in ([TB|), has non vanishing spin S = ^. 
In this case a further operator is relevant to generate all the manifold of degenerate states 



since by applying (m ■ S) to an element \S,m'^,p > of the tower ( |TB| ) one gets a different 
state with the same spin (due to the rotationally invariant expression [rh ■ S), commuting 
with the total spin) but momentum p + Q: 

\S,m'^ ,p + Q >(x {rh ■ S)\S,m'^ ,p > . (14) 

In fact, by using (p!TD, the operator {rn- S) commutes with m"*" and obviously annihilates any 
singlet state. Hence the expression (|14D defines non vanishing states only in the half-odd 
integer spin case where the lowest spin state of the tower has S = ^. The new relation (p!4D 
cannot be iterated as the previous one (|T^), because the square of the operator {S ■ m) when 
applied to the generator state \S = \, rn?,p > behaves as a constant equal to + o( j). In 



fact (m ■ S) 
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+ -^2- and S^S^ + = ^Sij on a spin-i state. 
The final diagram of the tower of degenerate states (or almost degenerate states at finite 
size) are shown in the following diagram: 



\S + l,m^,p + Q > # \S + l,m^,p> 
{S ■ m) 

\S,rn?,p> ^ \S,rn?,p + Q> 
{S ■ m) 



(15) 



We expect therefore that for odd number of holes and fixed spin (and maximum azimuthal 
spin) all the spectrum is characterized by couples of eigenstates with momenta p and p + Q 
related by the approximate expression (|14|). 

Indeed for the single hole case it can be rigorously shown (see App. |Vllip that the lowest 



eigenstates with momenta p and p + Q differ in energy by terms less than ~ j. Moreover 



the relation (|Tj) between such a couple of states is exact within the assumption that in any 
subspace with given spin and momentum the gap to the first excited state scales as ~ j ( 
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one has to excite at least one k ^ 0,Q magnon in this case), i.e. much larger than the energy 
accuracy ~ of the state ([1^). In the next section ( [V^ it will be shown that Eq. (|14D is 
indeed consistent with the exact solution obtained for momenta p = and p = Q at the 
supersymmetric point J = 2t. ||T3| 



The general character of the spectrum in a quantum antiferromagnet does not depend 
upon doping if long range magnetic order exists with finite momentum Q. However the basic 
relation ([14D is expected to be only approximate at finite doping, because charge excitations 
with energy cost ~ j are known to exist (invalidating the assumption of a gap scaling as 
j), as it is also clear from the weak coupling theory in the Hubbard model. 

It is a luck that for a single hole, charge excitations are forbidden by the requirement 



of fixed total momentum, and relation (|1^) is asymptotically exact for the lowest possible 
energy in each subspace with given momentum. Lowest energy eigenstates with momenta 
differing by Q are characterized by quite different eigenf unctions satisfying relation (p!4D. In 
two dimension momenta p and p+Q are equivalent by spatial symmetry on the surface of the 
magnetic Brillouin zone. Thus one expects some discontinuity or at least some singularity 
as the momentum of the hole crosses the magnetic Brillouin zone. This should hold at least 
for physical quantities, like the quasiparticle weight (see Sec. |V|), explicitly depending on 
the momentum dependent lowest energy state {ipp >■ 

IV. EXACT INTEGRATION OF THE SINGLE HOLE CHARGE 

The hamiltonian (^ is translation invariant and the most general one-hole state with 
total lattice momentum —p ( hole momentum p ) can be written: 

\^P>= -^T.^''"'fhTR\S> (16) 

where 15* > is a pure spin state without holes, i.e. fi\S >= 0. 

The state {ipp > in (|1^) represents the most general one- hole state {Nh = 1) with given 
hole momentum p. The action of any translation invariant operator O over the state {ipp > 
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does not change the momentum of the resulting state 0\il)p > and is therefore equivalent to 
the action of an effective spin operator ©'^■^■^ acting on 15 > , defined by: 

d\i;, >=^Y.^''''fRTRO^"\S > (17) 
For instance consider the kinetic term of the hamiltonian along one particular direction 

-^M = H fR'+T^fR'XR',R'+T^- (18) 
R' 

Using that the state 15 > and all the translated ones T^IS > contain no holes fRif\{TR\S > = 
^r,r'Tr\S >, K^\ipp > is easily computed: 

K,\ijp >= ^ E fR+r,XR,R+r,e^'"'Tn\S > 

V-tv R 

Then changing R + R in the dummy summation and using the translation operator 

rule TR_r^ = TnT.r^: 

y^R y^R 

where in the latter equality T_jiXr-t^,rTji = X-t^,o comes directly from (^. Finally, 
consistent with Eq. (p!7|) , it follows that: 

= Xo,-r,e-''^^T.^, (19) 
Analogously it is a simple algebra to show that the effective spin hamiltonian reads: 

H;" = E XCrM'^^T,. - J/^) + Hsw (20) 
where Hsw is the translation invariant Heisenberg hamiltonian 

Hsw = 7 E XR,R+Tf,- 

The total spin 5**°* = J2 ^Rof the effective spin hamiltonian remains unchanged with respect 

R 

to the old definition (H) and measures the total spin in the whole lattice including the origin 
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site O, whereas the pseudo-spin L, the physical spin S and the staggered magnetization rh 
turn in: 



L-ff _ 


So 






, Qtot Q 

O — Do 


(21) 


m^ff - 


So 

^ m — 

ij 


(22) 



The pseudo-spin and the physical spin operators commute with the effective spin hamil- 
tonian (pO|), as it is easy to check. Thus there is a one to one correspondence of any eigenstate 
of the effective spin hamiltonian with any single-hole eigenstate with given momentum of the 
extended hamiltonian (^. The identification of the true eigenstates of the t — J model is in 
this case trivial, because the total pseudo-spin L is fixed to | and the extended hamiltonian 
does not have eigenstates with unphysical pseudo-spin L ^ This property is valid only 
for the one hole it is discussed in [HI. 

In the following it will be established a correspondence between the eigenstates of the well 
known Heisenberg hamiltonian and the eigenstates of H'^^f . To this purpose it is convenient 
to use the total spin 5"*°* and its z— component S**"* as good quantum numbers for the 
eigenstates \S^°\ 5"*°* > of (^. The total spin of the effective hamiltonian S*°^ = L^^f + S^^^f 
corresponds to the sum of the physical spin with the pseudospin. It commutes with the 
hamiltonian but does not commute with each component of the pseudospin So- 

In order to have a definite total spin S"*"*, each eigenstate with physical spin S (with all 
the degenerate 2S + 1 components) has to combine with the two values of the spin at the 
origin, yielding eigenstates, with total spin S*"* = S + ^ and with S^°^ = 5 — |, by the well 
known addition relations of angular momenta. After projecting each eigenstate 15"*°*, Sl°^ > 
onto the ones with definite Sq = |, i.e \So >= (| + 5^)1 5*°*, 5^*°* > the physical spin of 
the corresponding hole eigenstates (|T6|) is 15**°* ± 1/2|. The physical spin is thus univocally 
determined to be S' = | in the singlet total spin 5**°* = subspace. It is possible to avoid 
any ambiguity by restricting all the following analysis to this subspace, which in turn is 
important for the analysis of photoemission experiments in stechiometric compounds. 
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In fact these experiments determine the imaginary part of the Greens function, that, at 
half filling, reads: 

G{p, t) = -i< H\cl^e-'^^ ^o)^Cp,.|i/ > (23) 

, i.e. it is obtained by creating a hole over the singlet antiferromagnetic state \H >. |T1 
Thus only the singlet subspace of the effective hamiltonian is relevant for the calculation of 
the Greens function in an antiferromagnetic insulator. 

In the singlet subspace we can use the total spin 5*°* to classify the spin of the elementary 
excitations, analogously to what was done in the ID Heisenberg model where the spinous 
have been found to carry spin i. [|15 



The hamiltonian ( pOD is exact , and the presence of the translation operator makes diffi- 
cult to use standard approaches as for the simpler Hsw- The recent ansaz proposed by 
Shraiman and Siggia corresponds to a variational semiclassical solution of the hamiltonian 
(pop, yielding for example the Neel state for Hsw- However in the one hole case the semi- 
classical solution cannot be controlled by the small parameter 1/s, important to derive the 
spin-wave limit for the Heisenberg model. |TB 



V. QUASIPARTICLE WEIGHT, GREENS FUNCTION AND CURRENT 

OPERATORS 

After the introduction of the effective spin hamiltonian (|20|) the Greens function (|23|) is 
easily expressed as an expectation value of a spin operator acting on the Heisenberg ground 

state \H >. 

The state \ipp >= Cp^a\H > is of the form (p^, if we choose 15" >= jS'/^ > with: 

\Sh >= n^,o\H > . (24) 

Due to the correspondence of eigenstates between -f/^'^'^ and H, we can expand 15* > in terms 
of eigenstates of H^ff and easily check that the propagation of \Sh > with the effective 
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Hamiltonian, \Sh >t= e p \Sh >, corresponds exactly to the propagation of i/jp with 
the exact t — J Hamiltonian and the Greens function immediately follows: 

G(p.() = -i<S«|e-'(*?''-«)'|S„> (25) 

Using that \Sh >= V^ni^o-\H >, that the commutator H^^^^Uo-^o 



vanishes and that G 



does not depend on a, we get, after Fourier transform G{p,uj) = J dt G{p,t) , 



G{p,u;) = l< H\ ^ \H > . (26) 
^ uj + 10 — Hp-' ■' 

The factor | is usually omitted in the literature of strong coupling theories like the t — J 
model, probably for estetic reasons of normalizations. Here we are interested to the actual 
photoemission delta weight and we are not allowed to use misleading normalizations. 

Once Hp-^f is diagonalized by eigenstates \i > with total vanishing spin S^°^ — S -\- Sq 
and energies Ei the Greens function is obtained by inserting this complete set of eigenstates 



in (H): 



G(p, u;) = W\<H\p>. \' l p (27) 
2 ^ uj + 10 - hi 

A general relation satisfied by the Greens function of a single hole in an antiferromagnet 
directly follows from the property that the lowest eigenstate \p + Q > oi H^^q can be 
written in term of the eigenstate \p > of H^-^-^ using relation (0) and (^) valid at low 
energy, i.e. \p + Q >iOC {S'^-^^ ■ rff^^)\p >j= —(So • frf^^)\p >i= —{Sq ■ m)\p >j where in the 
latter equalities we have used that \p >i are singlet states and we have neglected the 0{j) 
difference between the physical staggered magnetization m'^^^ ( P^ and the one acting over 
all the sites m = m'^-^^ + A simple normalization is then possible using Eq. (^Tp, and by 
neglecting 0(-^) contributions one obtains consistently: 

\p + Q >i= -{So-m)\p>i . (28) 
m 

for the lowest energy state in each fixed momentum sector. 
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The imaginary part of the Greens function is experimentally accessible by angle resolved 
photoemission experiments (ARPES) and we will refer to it as the spectral function A{uj). 
In the imaginary part of (^) some of the 6 function weights may remain finite for L —>■ oo 
and define the so called quasiparticle weight Zp, which is usually a single peak located at 
the bottom of the finite size spectrum |[T^. All the other part of the spectrum merge in a 
continuum of states for L — oo leading to an incoherent spectral function. 

Using (^TD the quasiparticle weight is given by: 

Zp = l\<H\S>p\' (29) 

where IS* >p is the lowest energy singlet state of the hamiltonian H^^^ . Notice that the 
quasiparticle weight in a strong coupling theory cannot exceed the value \ since it is obviously 
bounded by the value of the momentum distribution Up =< H\c^pCp\H >= |. Notice that 
using (p8D Zp^Q can be expressed in the following form: 

2 



<H\{m-So)\S>p\' (30) 



This is an asymptotically exact relation for L oo, as long as m > 0, i.e. within the 
assumption of long range magnetic order. 

Another important quantity to study is the current operator, which is useful when we 
calculate the transport properties. On a discrete lattice, the current operator corresponding 
to a uniform field is defined by ||18|: 



ietY^C^R^CR+r^^ + h.C. 



(31) 



In terms of spinless fermion hole operators fi and spin interchange ones Xij the total current 
is given by: 



J, 



iet fRfR+r^XR,R+r^ + h.c. 



RiT 



(32) 



We note that = —ie{K^ — K^^) where is defined in ( pTSD ,thus when applied to 
a state of the form (|TB|) the effective current operator , acting only on a spin wavefunction, 
can be written as 
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Jf = [ietxor/''^^T,, + h.c] (33) 



Following |]T9l the real part of the conductivity in a uniform field is characterized by a 



5 function proportional to the hole kinetic energy per site and a paramagnetic contribution 
given by the Kubo formula: 

2 ^M + A^(^) 

where = tXo,T^e*^'^'' is the kinetic energy in the direction fi and A is given by: 

ImA^(u;) = XI I < P\Jf^\P >j I'^^i^ - ^j) (35) 
j 

where the sum over j indicates the eigenstates with he same momentum p of the single hole 
ground state, uij being the corresponding energy excitations referred to the ground state 
energy. 

Of course a factor have been dropped out in the previous expressions for the conductiv- 
ity because there is only a single charge allowed to move. At small doping the conductivity 
will be proportional to the number of carriers times the conductivity of the single hole, i.e. 
will be finite in the thermodynamic limit. 



VI. SOME EXACT EIGENSTATES FOR THE SUPERSYMMETRIC POINT 

At the supersymmetric point for J = 2t the kinetic part proportional to t is exactly 
canceled by the magnetic bonds around the origin in (PPI). In fact any eigenstate |g > of 
the Heisenberg hamiltonian Hsw with total momentum g = p is an exact eigenstate of (^Of ) 
with the same hole momentum p and energy Ep. In fact by assumption Hsw\p >= Ep\p > 
and Tt-^Ip >= e~^^'^^'-\p >, and immediately follows that: 

H;ff\p >= Ep\p > +tY^xo,rM'^'T^. - l)b >= Ep\p > (36) 

The collection of all possible one hole states with arbitrary hole momentum has clearly 
an Hilbert space dimension ~ L times the dimension of the Heisenberg model. Thus the 

17 



Heisenberg eigenstates with all possible momenta q are not a complete set for the one hole 
Hilbert space. It is remarkable however that for J = 2t a considerable fraction (~ -^) of 
all the eigenstates is exactly known and indeed coincide with the ones of Hsw- I have 
verified numerically up to 26 site 2D-lattice that the lowest-energy one-hole eigenstates with 
momentum p = or p = Q = (vr, vr) are of the previous type, i.e. particular eigenstates 
of the Heisenberg model. For p = in any spatial dimension the singlet spin state IS* > 
characterizing the single hole eigenstate (|l^) coincides with the true ground state of the 
Heisenberg hamiltonian \H >^ the quasiparticle weight is exactly one and exhausts all the 
spectral weight, Ap=o = \^{!^ ~ Ep=o). 

An interesting feature instead is when the momentum of the hole coincides with the an- 
tiferromagnetic wavevector Q = (vr, vr, ■ ■ ■). In this case the singlet state \S > characterizing 
the lowest possible one hole eigenstate with momentum Q can be obtained using the first 
excitation of the Heisenberg model, which has momentum Q and is a triplet: 

\H, a > for (T = -1, 0, 1 
\H,1>= ^S^\H,0> \H,-1>= ^S-\H,0> 

with energy Ep=Q + + o{\), where xst is the static spin susceptibility of the Heisenberg 
model. The above states can be combined with the operator Sq which commutes with the 
hamiltonian, leading to nine different spin eigenstates Sq\H(j > for j,a = 1,2,3 with the 
same energy and with spin components on the S = 0,1,2 subspaces. By a proper linear 
combination of these states it is then possible to select an exact singlet eigenstate in the 
following way: 

1^ >Q= a (^S^S-\H,0> -SoS+\H,0> +2S^\H,0>) 
where the normalization constant is easily computed, yielding: 

After a little algebra, using that both \H > and cr > have definite momenta, the 

18 



overlap of this state with the Heisenberg ground state is: 

< H\S >Q= 6a < Hlm'glH, > 

, leading, by means of the Schwartz inequality, to: 

Zq = 1< H\S >l< 18«2 < H\{m^Qf\H >= 2m' (37) 

where the latter equality follow from \H > being a singlet, yielding < (m^)^ >= ^m'. 

The above inequality (^7|) is an upper bound for Zq, but represents an exact equality 
in the infinite size limit. In fact after a little algebra the state (^) with \i >= \H > 



saturates the bound (|37D , and represents the true eigenstate with momentum Q within the 
only assumption that the finite size gap in each sector of definite total spin and momentum 
scales as y (see App. |VlIl B|) . This basic assumption has been verified numerically as it is 



shown in Fig. 1 for t = and j = 2 for various momenta. In particular in the static limit 
t = the spin wave prediction of this gap always underestimates the value of the true gap 
computed by exact diagonalization on small lattices. In spin wave theory this gap scales 
as J since the perturbation induced by the hole affects only by small shift the bare spin 
wave energy excitation e^. The dispersion depends linearly with momentum ~ c\k\, 
where c is the spin- wave velocity, and the lowest gap is given by c^^stt ~ A plot of 
the quasiparticle weight for momentum Q = (vr, tt) is shown in Fig. 2 indicating that the 
finite size estimate of this quantity is strongly size dependent and without an exact result 
Zq = 2m' it would be difficult to decide whether Zq remains finite in the infinite size limit. 
This is in general the case for J ^ 2t, where the situation is still unclear and controversial 
so far. 

The exact determination of the quasiparticle weight, obtained for particular wavevectors 
and J = 2t has evidenced a more general property of the Greens function in an insulator 
with long range magnetic order. In fact hole momenta differing by the antiferromagnetic 
wavevector Q are in general characterized by the same energy in the infinite size limit, but 
with a quite different quasiparticle weight. The two one hole eigenstates with momenta p and 
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p + Q inside or outside the magnetic Brillouin zone respectively, have in fact a substantially 
different overlap with the Heisenberg ground state, because they are essentially derived from 
two orthogonal states: the Heisenberg ground state, and the corresponding lowest triplet 
excitation, respectively. This property is clearly a general one, valid for all momenta: In 
fact the relation between the lowest singlet eigenstates of H!^^^ can be written, analogously 
to Eq. (HD, : 

2 ^ 

IS" >p+Q= —{So ■ rn)\S >p 

At the boundary of the magnetic Brillouin zone, momenta differing by Q are equivalent by 
spatial symmetries. Then arbitrary close to this surface there should be at least a singularity 
in the quasiparticle weight, because the lowest energy state change dramatically even with an 
arbitrary small variation of momentum. It is reasonable to expect a jump of the quasiparticle 
weight that according to ( |37D is given approximately by: 

£(|±P^(.„)^^0.3T 
where we have used that in the Heisenberg antiferromagnet the numerical value for the 2D 



order parameter is m ~ 0.305. 



VII. CONCLUSIONS 

In this work I have discussed a general property of quantum antiferromagnets that can be 
detected by well resolved photoemission experiments that are determined by the low energy 
dynamic of a single hole in the undoped material. 

Contrary to the existing folklore considering the single hole problem as an old, boring 
and solved issue it is found here that at least an effect has been overlooked by the 

previous literature, effect that has been recently observed in the photoemission experiments 



by Wells et al. and were independently predicted in 

Such experiments have evidenced the surprising effect that the photoemission spectra of 
a quantum antiferromagnet apparently show up the presence of a Fermi surface as in the 
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corresponding metal at finite doping. 

In the present work this fact is a consequence of a zero-energy magnon excitation carrying 
the antiferromagnetic-wavevector momentum and having an infinite lifetime in an antifer- 
romagnet. This magnon is present in the ground state of the hole for momenta outside 
the Brillouin zone , while is absent in the other momentum region. The matrix elements 
entering in the spectral weight measured experimentally is jumping discontinuously along 
the mentioned surface separating the two momentum regions. This effect is just indicating 
the presence of this anomalous excitation. 

This kind of excitation is washed out in any mean field treatment since as discussed in 
the introduction, in such a case the direction of the order parameter is fixed and the spin is 
no more a measurable quantity. 
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VIII. SOME EXACT RESULTS FOR S =1 



Consider the effective hamiltonian (|20|) for the single hole problem: 

Hp = Hl + H' (38) 
Hl = tY.e'^^-Xo,rTr, (39) 

^' = ^ E {XR.,R,-\) (40) 

jjj represents the translation invariant Heisenberg hamiltonian without all the bonds con- 
necting the origin of coordinates O. The single hole hamiltonian commutes with the total 
spin: 
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S'"' = Y.Sr (41) 

R 

and the spin at the origin Sq, whereas the operator Sq = e~'''^^Sji defines long range 
order on a state \ip > if < iplS^Q ■ Sq\iP >^ im?L for L — > oo, where L is the number of 
sites, and m > is the value of the order parameter. 

We want to prove the following theorem: 
Theorem. Given a single hole eigenstate \p > of the hamiltonian Hp with 5**°* = 0, and energy 
Ep, assuming that long range order exists in the given state for momenta Q = (vr,??, . . .), 
corresponding to antiferromagnetic long range order, then the triplet states: 

|^^>=4b> j = 1,2,3 (42) 

define a state of the hamiltonian -f/p+g with energy expectation value 

^ <iIj^\iIj^> ^ L ^ ^ 

where the constant A is given by: 

-4Ep + 2<p\ETr,e'P^^\p> 
A = ^ . 4^ (44) 

where —Cp < is the ground state energy per site of the Heisenberg model obviously 
independent of p. The latter limit is easily obtained, from the definition (|4^ ) since T^-^ 
has all eigenvalues bounded by one. 

Proof. Since |p > is a singlet < ip^ip^ >= ^Lm^ is independent of j, as well as Ej. Then 
we consider the following operator: 



^ = E = ^ E [ShiHp^QSh - S^Hp) + {S^qHp^q - HpS^Q)Sl,] (45) 
By the exact relation: 



ShHp+QSi, = \{HpS^qS^q + Si,Si,Hp) + (46) 



using that , by assumption, Hp\p >= Ep\p > it easily follows that: 
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< p\F\p > 



Moreover from Tr^Sq = —SqTr^ and: 



o3 



<p\Sq ■ S-q\p > 



( q q qj qj \ 



yielding for S* = 1/2: 



8 

[^Q^ XR^R-i^Siq = -j{xr„Rj - 1/2), 



and the following expression for F holds: 



3 "^M 



(47) 



(48) 



4 9f 



Finally using the above two relations and the fact that Ej is independent of j the state- 
ment is easily proven (|43|). 



A. Consequences and remarks 



Suppose that the state with momentum p is the lowest energy state within the restriction 
of momenta p and p + Q in the Brillouin zone, and suppose that this lowest energy Ep is 
obtained when the total spin is minimum S* = 0, it will be shown that normalized singlet 



state defined in terms of m = -^Sq: 



\p + Q >= aL^So ■Th)\p> 



(49) 



has an energy expectation value Ep+g =< p + Q\Hp^Q\p + Q > arbitrary close to Ep and 
yielding the following bounds for the lowest energy Ep^q with momentum p + Q in the 
singlet subspace: 

3A 



Er, < Er,+o < E„ + — — 



(50) 



Proof. From: 
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{So ■ mf = + i(5o ■ m) - ^(5o ■ 5*"*) (51) 

the normalization constant is given by = ^ + 0(-^) (the expectation value of {Sq ■ m) 
on the singlet state \p > can be bounded using the Schwartz inequality for fixed component 
j: < p\{So ■ rn)\p >— 3 < p\Sjjm^\p >< ^m) Consider now the normalized state: 

\j >^ V3aLSi>m^\p> (52) 
This state has non vanishing overlap with the corresponding singlet one \p + Q >'■ 
«o = I < Jb + ^3 > r > (Re <j\p + Q>f^^ + O(^). 
Since the hamiltonian commutes with Sq < j\H\j >— — 7§r^ — which corresponds 

<P\SqSq\p> 

exactly to Ej in the previous theorem. By applying the theorem 

<j\H\j>^Ep + A/L 

On the other hand the state \j > can be written as \j >— ao\p + Q > > where the 
state i/j' has no component in the singlet and < ip'lip' >= 1 — a^. By hypothesis this state 
has an energy expectation value — <^^|^'?> — higher than Ep, the minimum possible energy 
between the p and p + Q subspaces. Thus using the conservation of the spin, it follows that: 

Ep + A/L =< ii\H\ii >> alEp+Q + (1 - al)Ep 

where -Ep+g =< p + Q\H\p + Q > is a variational estimate of -Ep+Q. Finally it the latter 
inequality gives: 

Ep < Ep^Q < Ep^Q < Ep + 3A/ L 

which concludes this proof. 



Remark. The theorem is more generally valid even when the lowest energy state is no more 
a singlet state provided the energy gain to the S* = 1 and S = 2 vanishes at least as 2;. In 
fact the trial states \ j > have component only in the S — 0, S — 1 and S — 2 sectors. 
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B. Exact variational state 



It has been shown in the previous section that the variational state \p -\- Q >, defined 
in the singlet subspace is arbitrarily close in energy to the exact lowest energy singlet state 
\ipp+Q > with momentum p + Q. It is reasonable to assume that the finite size gap in each 
subspace with definite spin and momentum is of the order ^ because it is determined by 
an excitation of at least one magnon (remind that the spin wave excitations are of order 
c\k\, where c is the spin wave velocity and the minimum allowed \k\ > is of order j). 
If the above hypothesis is correct as it can be easily verified numerically (see Fig. 1) it is 
possible to show rigorously that the variational state \p + Q > is arbitrarily close to the 
exact eigenstate ipp+Q- 

In fact suppose \p + Q >— ao|'0p+Q > +W where by definition > is orthogonal to 
the lowest state with definite spin and momentum, and thus satisfying by assumption: 

< 'il)'\H\'il)' > 



Then it easily follows that: 



<p + Q\H\p + Q>= alEp+Q + (1 - al){Ep+Q + B/l) = E^+q + (1 - al)B/l. 



The previous relation is compatible with the previously stated theorem only if 

1-al =<V;'|V^'><||-0 



which proves the statement of this section for d > 1, as L = Z 



d 



25 



REFERENCES 

[1] B.O. Wells et al. Phys. Rev. Lett. 74, 964 (1995); Z.-X. Shen, W. E. Spicer, D. M. King, 
D. S. Dessau, and B. O. Wells, Science 267, 343 (1995). 

[2] P. Aebi et al. Phys. Rev. Lett. 72, 2757 (1994). 

[3] S. Chakravarty et al. Phys. Rev. Lett. 60, 1057 (1988). 

[4] A. Parola and S. Sorella, Phys. Rev. B 45, 13156 (1992); S. Sorella and A. Parola, J. 
Phys. Cond. Mat. 4 3589 (1992). 

[5] C.L. Kane, P.A. Lee and N. Read, Phys. Rev. B 39, 6880 (1989), S. Schmitt-Rink, C. 
M. Varma, and A. E. Ruckenstein, Phys. Rev. Lett. 60, 2793 (1988). 

[6] G. Martinez and P. Horsch, Phys. Rev. B 44, 317 (1991). 

[7] Z. Lin and E. Manusakis Phys. Rev. B 44, 2414 (1991). 1 

[8] A. Angelucci, S. Sorella and D. Poilblanc, Phys. Lett. A 198 , 145 (1995). 

[9] B.I. Shraiman, E.D. Siggia Phys. Rev. B 42, 2485 (1990). 

[10] see e.g. F. Strocchi "Elements of Quantum mechanics of infinite systems" World Scien- 
tific, Singapore (1985). 

[11] H. Neuberger, T. Ziman, Phys. Rev. B 39, 2608 (1989). 

[12] B. Bernu, C. Lhuillier, and L. Pierre Phys. Rev. Lett. 69, 2590 (1992). 

[13] P. A. Bares, G. Blatter and M. Ogata Phys. Rev. B 44 , 130 (1991). 

[14] E. Lieb and D. Mattis J. of Math. Phys. 3, 749, (1962). 

[15] L. Faddeev, L. Takhtajan, Phys. Lett. A 85, 375 (1981). 

[16] P.W. Anderson, Science 235, 1196 (1987). 

[17] Q. F. Zhong and S. Sorella Phys. Rev. B 51 16135 (1995). 

26 



[18] D. J. Scalapino S. R. White and S. Zhang Phys. Rev. B 47, 7995 (1993). 
[19] F.C. Zhang and T.M. Rice, Phys. Rev. B 37, 3759 (1988). 
[20] Q.F. Zhong, S. Sorella, Europys. Lett. 21, 629 (1993). 
[21] S. Sorella, J. Low Temperature Phys. 99, 371, (1995). 



27 



FIGURES 

FIG. 1. Plot of the finite size gap to the first singlet excitation as a function of the inverce lattice 
side length in 2D for (a) the static limit t = , full dots arc obtained by numerical diagonalizations 
and triangles represent the spin- wave predictions, (b) J = 2t and momentum p = (0,0). In this 
case also the gap to the lowest S = 3/2 excitation of the Heisenberg model scaling to zero as 
is shown for comparison (squares), while the triangles represent the corresponding gap for the 
Heisenberg lattice (without the hole) in spin-wave approximation, (c) J = 2t and momentum 
p = (7r/2,7r/2) corresponding to the single hole ground state. The triangles are the same points as 
in (b). (d) J = 2t and momentum (7r,7r). The triangles are the same points as in (b). 

FIG. 2. Plot of the quasiparticle weight for momentum p = (vr, vr) as a function of the lattice 
size for different values of j. The star indicate the expected asymptotic L —>■ oo value (see text) 
for J = 2t. 
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